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ABSTRACT 

This  work  deals  with  the  correlation  functions  and  their 
methods  of  obtention.   A  method  without  a  "pure"  time  delay 
is  considered  which  uses  Laguerre  functions  type  filters. 
The  possibilities  of  utilization  of  the  nonsymmetric  Laguerre 
type  filter  are  analyzed.   A  very  inexpensive  set  of  circuits 
which  are  useful  to  obtain  the  mentioned  filters  for  low 
frequencies  signals  are  presented,  together  with  their 
computer  analysis  and  laboratory  realizations. 
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I.   INTRODUCTION 

Correlation  analysis  has  become  the  most  powerful  tool 
to  separate  the  signal  from  the  noise,  especially  in  those 
applications  where  the  signals  are  hidden  beneath  a  blanket 
of  other  signals  and  extraneous  noise. 

The  development  of  instruments  which  can  continuously 
sample  even  the  most  noisy  signals  and  compute  the  correla- 
tion function  simultaneously,  have  extended  the  use  of 
correlation  analysis  into  radio  astronomy,  fluid  and  solid 
state  physics,  neurology,  seismology,  geophysical  exploration 
and  other  areas  [Ref.  3]. 

Radar  and  sonar,  of  primary  importance  in  military  appli- 
cations, use  correlation  detectors  because  it  has  been  shown, 
for  the  case  of  the  binary  decision  problem,  that  no  matter 
what  the  philosophy  employed  for  the  design  of  a  receiver 
is,  the  best  detector  is  always  a  crossconrelator  detector 
[Ref.  5]. 

Theoretically,  an  infinite  signal  to  ncoise  ratio  can  be 
obtained  in  the  detection  of  periodic  signals  in  noise 
through  the  correlation  technique,  but  this  would  require 
an  infinite  time  of  observation.   In  practice,  by  using 
crosscorrelation  practical  measurements  can  be  made  in 
reasonable  time. 

Initially,  theoretical  considerations  on  correlation 
functions,  their  properties  and  methods  to  obtain  them  using 


pure  time  delays  will  be  presented.   Next,  this  thesis  will 
present  a  method  which  does  not  use  a  pure  time  delay  and 
yields  results  comparable  to  classical  methods. 

The  mathematical  formulation  of  the  method  and  the  cir- 
cuits that  can  be  used  are  studied  and  compared  with  the 
standard  methods.   Finally  an  approximation  to  the  necessary 
type  of  circuits  was  analyzed  by  digital  computer  methods 
and  checked  in  the  laboratory  using  actual  components. 


II.   CORRELATION  FUNCTIONS 

It  is  assumed  that  some  physical  process  produces  the 
time  functions  x..(t),  Xp(t),  ....  x^Ct)  simultaneously. 
It  is  assumed  further  that  the  physical  process  is  stationary, 
in  other  words  it  is  not  changing  over  time.   Also  it  is 
assumed  that  the  time  functions  are  not  zero  and  do  not  have 
a  d.c.  component  or  they  have  zero  mean  if  the  time  functions 
correspond  to  a  stochastic  process  [Ref.  3]. 

A.   AUTOCORRELATION  FUNCTION 

The  time  average  autocorrelation  function  for  determinis- 
tic or  known  signals  is  defined  for  two  different  classes  of 
signals,  signals  of  finite  average  power  (case  II)  and 
signals  of  finite  energy  or  pulse  signals  (easel). 

For  signals  of  case  II  the  definition  is: 


R- 


(t)  =  lim  ^t   /   x,(t)  xn(t  +  t)  dt        (1) 
T-m»  CL    -T 


and  for  signals  of  case  I: 


(£(t)  =    /"xjU)  xx(t  +  t)  dt  (2) 

—00 


If  the  time  function  under  consideration  is  a  stochastic 
process  the  autocorrelation  function  is  defined  as  follows: 


R(t)  =  E[x1(t)  xx(t  +  t)]  (3) 


If  the  stochastic  process  is  ergodic 

(R(t)  =  R(t)  (4) 

According  to  the  definitions  given  above,  Figure  1  shows 
an  elementary  device  to  obtain  the  autocorrelation  function. 
The  process  implies  multiplication  of  the  signal  by  its 
delayed  replica  and  averaging.   If  the  instantaneous  values 
of  x..(t)  and  x,  (t  +  t)  are  multiplied  and  the  product,  averaged 
over  a  sufficiently  long  time,  the  result  is  the  autocorre- 
lation function  if  some  means  to  vary  t  have  been  added  to 
the  circuit. 

1.   Properties  of  the  Autocorrelation  CRef .  3] 

The  autocorrelation  function  has  the  following 
properties: 

a.  Maximum. 

The  autocorrelation  function  is  a  maximum  at 
t  =  0,  mathematically: 

|(R.(x)|  <&(0)    Schwartz  Inequality         (5) 

b.  Value  at  t  =  0 

The  value  of  the  autocorrelation  at  x  =  0  is 

connected  to  the  total  power  of  the  signal.   If  x..(t)  is  a 

2 
voltage,  then  for  t  =  0,  x..(t)   Is  the  average  power  of  the 

signal  (as  measured  on  a  1  ohm  resistor).   Hence\JC(0)  is 

the  average  power  of  signals  of  case  II  andvK(O)  is  the 

energy  In  the  pulse  signals  of  case  I. 
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FIGURE  1:   Typical  Correlation  Functions. 
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c.  Function 

The  autocorrelation  function  is  a  function  of 
the  form  of  which  is  characteristic  of  the  original 
gnal  x,  (t) . 

d.  Type  of  Function 

The  value  of  the  function  for  negative  values  of 
Is  identical  to  that  for  the  same  positive  values.   This 
ans  that  the  autocorrelation  function  is  an  even  function 
x: 


(R.(t)  =Uu-t)  (6) 


e.  Averaging  Time  \ 

If  the  averaging  time  is  long  compared  with  the 
ciprocal  of  the  lowest  frequency  in  the  original  signal, 
en  repeated  measurements  of  the  product  for  a  given  x  will 
eld  values  very  close  to  one  another. 

f.  Periodic  Functions 

If  x, (t)  contains  any  periodic  frequency  compon- 
ts,  thenvkld)  will  contain  each  of  these  frequency  compon- 
ts.   So  for  certain  time  f unctions, \K,(x)  is  known  if  the 
iginal  signal  x, (t)  is  given. 

g.  Physical  Processes 

Any  physical  realizable  process  that  produces 
time  function  like  x-.(t)  is  such  that  the  value  of  the 
notion  at  time  t  becomes  more  independent  of  the  value  at 
+  t  as  x  gets  larger.   Thermal  noise,  and  even  the  quantum 
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mechanical  uncertainty  principle,  introduce  randomness  that 
cause  a  gradual  loss  of  coherence  between  the  signal  and 
its  delayed  replica.   This  is  true  even  for  the  oscillators 
in  the  most  stable  atomic  clock. 

Hence  for  a  signal  arising  from  a  real  process 
the  autocorrelation  function  approaches  zero  as  t  becomes 
sufficiently  large. 

h.   Coherence  Time 

The  value  of  t  that  causes  a  significant  reduc- 
tion in  the  autocorrelation  function  is  a  measure  of  the 
coherence  time  of  the  original  signal  x,  (t). 
i.   The  Frequency  Domain 

A  less  obvious  property  of  the  autocorrelation 
function  is  indicated  by: 


(R.(t)<=»4>(u))  (7) 


Equation  (7)  means  that  the  autocorrelation 
function  and  the  power  spectral  density  are  Fourier  Transform 
pairs.   Mathematically: 

♦  (u)  =   /<R(T)e~j(0T  dx  (8) 

—  00 

(R.(x)    =  t£-       /°°<l>U)eJurr   du  (9) 

— oo 

The  power  spectral  density  for  deterministic 
case  II  signals  is: 
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♦  (a))  =  lim  2Y  x1T(w)  X^U)  (10) 


T-*» 


lim  2^r  I x1t(cd)  | 
T- 


where  the  subindex  T  indicates  that  the  function  x,  (t)  is 
truncated  and  hence  observed  through  a  window  of  observation 
of  duration  T. 

The  power  spectral  density  for  deterministic 
case  I  signals  is: 


<f>U)  =  X1(u)  X1*(u)  =  |X1(co)|2  (11) 


The  same  property  holds  for  stochastic  processes, 
in  which  case: 


and 


R.(t)^(m)  (12) 


♦  (a))  =  lim  ^r  E{  |X  (w)  |2}  (13) 

T-«» 


j .   White  Noise 

The  autocorrelation  function  for  very  wideband, 

uniform  (white)  noise  with  r.m.s.  value  E   is  an  impulse 

2 
function  at  t  =  0  with  amplitude  E  .   This  means  that  one 

characteristic  of  such  kinds  of  noise  is  that  its  instantaneous 

value  is  completely  independent  of  the  value  at  any  other 

instant  and  that  the  coherence  time  of  the  process  producing 

the  noise  is  practically  zero  or  very  short. 
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B.   CROSSCORRELATTON  FUNCTION 

While  the  autocorrelation  function  of  a  signal  is  equi- 
valent to  the  traditional  technique  of  power  density  spectrum 
analysis,  there  is  no  classical  analogy  for  crosscorrelation 
analysis.   Crosscorrelation  is  concerned  with  the  relation- 
ship between  two  different  signals  that  arise  in  some  common 
process  [Ref.  3]» 

The  mathematical  expression  of  crosscorrelation  is: 

(K-12(t)  =  lim  ^r   /   x,(t)  xg(t  +  t)  dt      (14) 
T-h»    — T 

1.   Properties  of  the  Crosscorrelation  [Ref.  3] 

The  properties  of  the  crosscorrelation  function  are, 
in  general,  quite  different  from  those  of  the  autocorrelation 
function. 

a.  Type  of  Function 

The  crosscorrelation  function  is  not  an  even 
function  of  t,  which  means: 

^12(l)  ^2("t)  (15) 

However,   (*.   (-t)  =^21(t)  (l6) 

a  relationship  that  has  practical  importance  in  obtaining 
\K2(x)  for  negative  delays. 

b.  Coherence 

The  crosscorrelation  function  approaches  zero 


15 


as  x  approaches  infinity  for  signals  that  arise  from  real 
physical  processes,  due  to  the  presence  of  noise  and  the 
uncertainty  principle . 

Furthermore  if  x..  (t)  and  x~(t)  arise  from  two 
completely  separate,  unrelated  processes 


(R12(x)  =  0  (17) 


c.   Frequency  Domain 

The  same  as  in  the  case  of  autocorrelation 
functions : 


a 


2(T)<^4>12U)  (18) 


In  this  case,  however,  the  physical  meaning  of 
$ ,p(aj)  is  not  so  clear.   It  can  be  called  the  spectrum  of 
crosscorrelation  of  the  time  functions  x,  (t)  and  x~(t) . 

Figure  2  shows  some  typical  examples  of  cross- 
correlation  between  functions. 

Crosscorrelation  analysis  provides  a  powerful 
analytical  tool.   The  ability  to  measure  the  degree  to  which 
the  signals  that  arise  from  a  common  physical  phenomenon 
resemble  each  other  as  a  function  of  the  delay  time  between 
them,  can  provide  a  much  deeper  insight  into  the  phenomenon 
being  studied  than  a  separate  analysis  of  the  properties 
of  either  signal  alone. 
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III.   APPLICATIONS  OF  CORRELATION 

A.   DETECTION  [REF.  4] 

Given  two  functions  consisting  of  a  signal  S(t)  and 
noise  N(t ) : 


Xl(t)  =  Sl(t)  +  Nl(t)  (19) 


X2(t)  =  S2(t)  +  N2(t)  (20) 


it  is  possible  to  crosscorrelate  them  with  the  following 
result: 


vK12(t)  -  11m  ^y   /   CS1(t)  +  N1(t)][S2(t+x)  +  N2(t+x)]  dt 


\s^    +«3  H  (t)  +\S2(X)  ««     (t> 


(21) 


All  the  terms  of  the  above  expression  will  be  zero  as 
X  approaches  infinity. 

Now  it  is  assumed  X, (t)  =  X«(t)  and  hence  the  four  terms 
of  the  crosscorrelation  become 


^11(T)  =^SS(T)  +^SN<T>  +<^NS(T)  +(RNN(T)       <22) 


(JL„CO  i: 


In  the  expression  thus  obtained  VK^^t  )  is  the  autocorre- 
lation function  of  the  signal  itself  and  is  a  function  of 
t  other  than  zero  even  for  large  values  of  x.   The  remaining 
terms  approach  zero  as  x  becomes  larger.   So  the  first  term 
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is  a  measure  of  the  signal  and  the  other  ones  are  a  measure 
of  noise . 

If  now  it  is  assumed  that  the  frequency  of  the  given  sig- 
nal is  known,  an  internal  generator  in  the  receiver  can  be 
used  to  generate  the  second  input  to  the  correlator  and  the 
two  signals  to  crosscorrelate  are: 


X1(t)  =  Sx(t)  +  N(t)  (23) 


X2(t)  =  S2(t)  (24) 


The  above  one  is  the  typical  case  of  a  radar  receiver, 
and  if  the  two  signals  are  crosscorrelated  the  result  is: 

Now  the  number  of  noise  terms  has  been  reduced  to  one 
and  it  is  reduced  to  zero  for  large  values  of  x.   The  pre- 
vious examples  show  how  correlation  techniques  can  be  used 
to  separate  desired  signals  from  the  noise. 

B.   ANALYSIS  OP  LINEAR  SYSTEMS  [REP.  3] 

It  can  be  shown  that  it  is  possible  to  obtain  a  unit 
impulse  response  of  a  linear  system  by  driving  it  with  broad- 
band-white-noise and  crosscorrelating  this  input  with  the 
system  output,  as  is  shown  in  Figure  3.   The  correlogram 
that  results  gives  the  same  information  about  the  system  as 
if  it  were  excited  with  an  approximation  to  a  delta  function 
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and  its  output  recorded  on  an  oscilloscope  or  x-y  recorder. 

This  technique  has  important  practical  applications 
because  of  its  immunity  to  internal  system  noise  and  because 
it  permits  to  obtain  the  impulse  response  of  a  linear  system 
in  the  presence  of  a  signal  in  the  linear  system. 

With  real  time  methods  of  correlation  it  is  possible 
using  the  method  described  above  to  keep  the  response  of 
critical  systems  under  constant  surveillance  and  to  make 
optimizing  adjustments. 
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IV.   METHODS  TO  OBTAIN  CORRELATION  FUNCTIONS 
There  exist  three  methods  to  obtain  correlation  functions: 

A.   METHOD  THROUGH  DEFINITION-TIME  DOMAIN 

The  easiest  way  to  obtain  a  correlation  seems  to  be  the 
process  indicated  in  its  mathematical  definition  (1)  and 
(1*0,  this  means:   a  delay  t,  is  given  to  one  of  the  signals, 
then  the  other  signal  and  the  delayed  one  are  multiplied  and 
finally  the  product  is  averaged  being  the  resultv(.(T, )  one 
point  of  the  correlation  function.   Varying  t  it  is  possible 
to  obtain  a  graph  representing  the  correlation  function. 

A  basic  block  diagram  to  perform  the  process  was  shown 
in  Figures  1  and  2.   The  essential  feature  of  such  devices 
is  the  temporal  delay  x  between  the  two  signals.   When  the 
signals  are  fluctuating  voltages  or  currents,  the  provision 
of  the  "pure"  delay  is  rather  a  problem  because  the  available 
delay  devices  produce  distortion  on  the  signals  for  long 
delay  times. 

Another  problem  of  this  method  is  that  it  is  essentially 
discrete,  and  the  number  of  points  of  the  correlation  function 
that  can  be  obtained  in  a  given  time  interval  depends  upon 
the  number  of  individual  delay  devices. 

Generally  to  avoid  distortion,  digital  delay  devices  are 
used  [Ref.  3]  and  hence  it  is  possible  to  obtain  a  number 
of  discrete  values  of  the  correlation  function  without 
distortion. 
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Although  this  method  is  not  cheap  it  is  widely  used  for 
receivers  which  use  correlator  detectors. 

B.  INDIRECT  METHOD-FREQUENCY  DOMAIN 

Equations  (8),  (9),  (10),  (11)  and  (13)  set  up  the 
relations  between  correlation  functions-time  domain,  and 
power  spectral  densities-frequency  domain. 

It  is  seen  from  the  mentioned  equations  that  correlation 
of  two  signals,  or  one  signal  with  a  shifted  version  of 
itself,  in  the  time  domain  is  equivalent  to  a  complex  con- 
Jugate  multiplication  of  their  linear  spectra  in  the  frequency 
domain . 

In  other  words,  the  autospectrum  and  the  autocorrelation 
are  related  through  the  Fourier  Transform,  the  correlation 
being  the  inverse  transform  of  the  power  spectrum.   The  same 
relation  applies  to  the  cross-spectrum  and  crosscorrelation. 

Although  this  method  is  the  faster  one  it  is  expensive 
because  to  obtain  correlations  means  have  to  be  provided 
to  obtain  the  power  spectrums  of  the  signals.   Hence  this 
method  is  typical  of  special  signal  processing  machines 
which  are  designed  to  observe  several  characteristics  of 
signals,  and  the  power  spectrums  are  not  only  used  to  obtain 
correlation  functions. 

C.  METHOD  WITHOUT  "PURE"  DELAY-TIME  DOMAIN 

This  method  has  the  advantage  that  a  "pure"  delay  is 
not  required.   The  correlation  is  obtained  as  a  summation  of 
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terms  and  the  method  is  based  upon  the  properties  of 
orthogonal  functions  [Ref.  2], 

This  is  the  only  method  essentially  analogous  and  as 
will  be  shown  in  the  following  sections  inexpensive 
correlators  can  be  built  for  some  special  kJLnd  of  signals 
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V.   CORRELATION  WITHOUT  A  "PURE"  TIME  DELAY 

A.   MATHEMATICAL  FORMULATION 
1.   Analysis  Problem 

It  Is  supposed  that  any  correlation  functionvK(x) , 
either  autocorrelation  or  crosscorrelation,  can  be  expanded 
in  a  series  of  orthogonal  functions  [Ref.  2],  thus 

(£(t)  =  E   an  6  Ct)[u)(t)]y  0  <  t  <  «        (26) 
n=0  n  n 


where  the  6  (t)  are  polynomials  which  form  an  othonormal  set 
with  respect  to  the  weight  function  oo(t)  in  the  range 
0  <_  t  <  «.   Then  the  coefficients,  a  ,  in  the  series  expansion 
can  be  found  in  the  usual  way  and  are  given  by 


CO—. 

an  =  /(K.(t)  9  (tHwCt)]1-7  dT  (27) 

n   0 


Now  the  system  of  Figure  4  is  considered.   The 
linear  network  has  an  impulse  response  h  (t) .   So  the  filter 
output  is : 


00 

xn(n;t)  =  /  x, (t-u)  h  (u)  du  (28) 

i        0 


and  the  multiplier  output  is 


xx(n;t)  x2(t)  =   /  Xl(t-u)  x2(t)  hn(u)  du   (29) 
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being  the  output  of  the  averaging  device 


x1(n;t)  x2(t)  =  /  xx(t-u)  x2(t)  ^(u)  du   (30) 


0 


/(R.(u)  h  (u)  du  (3D 


0        n 


Now  comparing  (3D  to  (27)  it  is  seen  that  the 
coefficient  a  is  identical  with  the  average  x, (n;t)  x~(t ) 
provided  that  the  filter  impulse  response  is  so  chosen  that 


hn(t)  =  9n(t)[o3(t)]1"Y  (32) 


The  filters  whoise  impulse  responses  satisfy  Equation 
(32)  are  called  "orthogonal  filters,"  and  iff  they  could  be 
built,  the  coefficients  a  ,  would  be  obtained  simply  by 
switching  in  each  of  the  filters  in  turn  and  determining  the 
corresponding  averages  [Ref.  2]. 

Naturally,  in  any  practical  system  only  a  finite 
number,  say  (N+l) ,  of  filters  would  be  used- 
2.   Synthesis  Problem 

Equation  (26)  can  be  written 

(R.  (t)  =  [cud)]2*""1   E   an  en(T)[a)(T)]1-Y      (33) 

0X1     II 


or 


(R.(t)  =  [Xt)]27"1  l     a  h  (t)   0  <  t<  oo     (3H) 


n=0  n  n 
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In  particular,  if  a  finite  number  of  terms  are 
taken,  say  N,  it  can  be  written  as  an  approximation  to  (33) 


CRM(t)  =  [a)(t)]2Y_1   Z   an  h  (t) 
N  n=0  n  n 


If  Y  is  taken  to  be  zero,  (35)  becomes 


(35) 


&„(t)  =  [wCt)]"1  E 

n=0 


N(t)  =  [co(t)]"x  L   an  hn(t)         (35a) 


and  if  Y  is  taken  to  be  1/2 


Q\*T(t)  =   z  a_  h  (t)  (35b) 


n=0  n  n 


Figure  5  shows  the  complete  synthesis  of  the  method. 

It  is  possible  to  use  half  the  number  of  linear  networks 

calculating  first  the  coefficients  due  to  the  particular 

signals  to  be  correlated,  and  then  after  introducing  the 

values  of  a  the  approximate  correlation  function  can  be 
n      ^ 

seen  on  the  screen  of  a  CRT  as  a  transient  following  the 
application  of  the  input  impulse  [Ref.  2]. 

A  convenient  system  results  if  in  fact  the  pulsing 
can  be  carried  out  repetitively  at  intervals  that  are  long 
compared  with  the  filter  time  constants,  and  yet  sufficiently 
frequently  to  give  a  steady  trace  on  the  CRT  screen.   This 
can  always  be  arranged  in  practice  [Ref.  2]. 

It  can  be  seen  from  Equation  (35a)  that  if  Y  is 
taken  to  be  zero,  the  output  of  the  circuit  in  Figure  5  has 
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to  be  multiplied  by  [w(t)]   to  obtain  the  correlation 
function. 


B.   TWO  REALIZABLE  SYSTEMS 

The  applicability  of  all  the  preceding  theory  depends 
on  realizing  a  set  of  filters  whose  impulse  responses 
satisfy  Equation  (32). 

The  polynomials  8  (t)  and  hence  the  impulse  responses 
h  (t)  depend  only  on  the  weight  function  aj(t)  and  the 
parameter  y . 

The  weight  function  u(t)  -  a  e    was  cihosen. 

1.   Case  of  y  =  0 

The  transfer  function  H(s)  of  a  single  R-C  cell 
(low  pass  filter)  is: 


s  +  - 

T 


being 


1    1 
a  '  7  -  RC 


HRC(s)  corresponds  in  the  time  domain  to  the  impulse 


response 


hRC(t)  =  a  e"at  (37) 


If  now  H(s)  of  a  single  C-R  cell  (high  pass  filter) 
is  considered: 

30 


hcr(s)  ■  -fr  ■  rhr  (38) 

s  +  - 


and  the  impulse  response  in  the  time  domain  is 


hCR(t)  =  -a  e"at  (39) 


Now  a  number  n  of  high  pass  filter  cells  are 
cascaded  and  between  each  cell  a  buffer  amplifier  is  con- 
nected.  Buffer  amplifiers  have  ideally  infinite  input 
impedance,  zero  output  impedance  and  unity  gain,  and  they 
isolate  the  cells. 

Hence  the  transfer  function  in  the  s  domain  for  the 
system  of  n  high  pass  filter  type  cells  is: 


H   (s)  =  [HrR(s)]n  =  — ^ (40) 

nCR        CR        (s  +  a)n 


If  a  single  low  pass  filter  type  cell  is  connected 
to  the  input  of  the  circuit  considered  above  via  a  buffer 
amplifier,  the  resultant  circuit  is  of  the  form  shown  in 
Figure  6  and  its  transfer  function  in  the  s  domain  is: 


H  (s)  -  H„p(s)  H  „_(s)  =   ? (41) 

n       RC     nCR      s  +  a   ,     ^n 


a  sn 


t~  xn+1 

(s  +  a) 


(s  +  a)J 

(42) 
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The  inverse  Laplace  transform  of  the  above  one  is: 


d    1  .n  -at  /in\ 

a  "TiT  HT  *  e  ("3) 

at 


which  can  be  written  as: 


hn.(t)  =  a  ."**  j  ^  fe^      (*« 


■  a  e~   L  (at) 
n 


where  a  e~   is  the  weight  function  and  L  Cat)  are  the 
Laguerre  polynomials  and  the  product  is  defined  as  Laguerre 
functions . 

The  given  expression  for  h  (t)  satisfies  Equation 
(32)  and  the  conditions  of  orthogonality  [Bef .  2]  and  hence 
the  circuit  of  Figure  6  can  be  used  to  buiid  the  correlator 
2.   Case  of  T  =  1/2 

Another  filter  based  on  the  same  weight  function  is 
the  lattice  structure  shown  in  Figure  7.   Im  this  case  the 
s  domain  transfer  function  is  given  by 

J  s  -  p  n 

Vs)  -  -hr   ( — l>  (il5) 

1  s  +  I   s  +  I 

and  the  impulse  response  is: 

-  °-t 

h  (t)  -  I  e  2  L  (at)  (46) 

n      2        n 
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which  is  identical  as  the  case  of  y   =  0  except  for  a 
constant  factor  [Ref.  2]. 

3.   Frequency  Domain  Considerations 

The  .nonsymme trie  filter,  case  of  Y  =  0,  was  used  with 
good  results  to  build  an  autocorrelator  [Ref.  2]  which  used 
N  =  10.   Hence  N+l  =  11  filters  were  used,  being  the  first 
one  a  low  pass  filter  (n  =  0)  and  the  rest  combinations  of 
a  low  pass  filter  and  a  different  number  of  high  pass  filters 
isolated  by  cathode  follower  amplifiers. 

We  saw  in  the  ideal  case,  which  corresponds  to 
perfect  isolation  buffer  amplifiers,  that  the  transfer 
function  in  the  frequency  domain  was  (4l): 

H(s)  "  r-r^  <rhr>n  <*« 

and  because  of  the  isolation  between  cells  it  was  taken  as 
the  product  of  the  transfer  function  of  a  low  pass  filter 
and  n  R-C  high  pass  filter  cells. 

For  the  Ju  domain  the  above  one  becomes: 

|H<J„)  |2  -  «2/"  5  +i  CH8) 

(a  +  a)  ) 

being  a  function  of  a  and  n. 

Figures  8,  9,  10  and  11  show  the  transfer  functions 
in  the  frequency  domain  of  a  R-C  low  pass  filter,  an  R-C 
high  pass  filter,  ten  high  pass  R-C  cells  cascaded  and 
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ideally  Isolated,  and  of  the  Ideal  filters  for  n  =  1  and 
n  ■  10,  respectively.   In  all  cases  a  =  10  was  taken. 

The  transfer  functions  were  obtained  with  the  aid  of 
the  CORNAP  circuit  analysis  program  [Ref.  1]  for  the  case  of 
individual  cells  and  for  the  case  of  multiple  cells  cascaded 
taking  the  product  of  the  individual  transfer  functions. 

From  the  curves  it  can  be  easily  seen  which  is  the 
frequency  range  for  the  incoming  signals  and  which  the  band- 
width of  the  system.   Also  if  the  incoming  signals  are  given 
it  is  possible  to  see  which  filters  need  postamplification. 

It  has  to  be  pointed  out  that  the  frequency  domain 
considerations  are  only  important  in  the  determination  of 
the  coefficients  (Figure  5)  only  case  in  which  one  of  the 
input  signals  is  convoluted  with  the  impulse  responses  of 
the  filters. 

As  the  system  is  linear,  by  varying  a  it  is  possible 
to  obtain  filters  for  different  ranges  of  frequencies. 

The  correlator  described  in  Ref.  2  used  a  =  0  type 
filters  to  autocorrelate  a  sine  wave  which  frequency  was 
152  Hz,  being  the  basic  time-constant  of  the  filter  a  =  843 
sec~"  .   The  results  were  practically  identical  to  those 
which  can  be  obtained  with  conventional  methods. 

C.   A  SET  OF  PRACTICAL  FILTERS 

The  buffer  amplifiers  needed  to  build  the  filters  have 
to  have  infinite  input  resistance,  zero  output  resistence 
and  unity  gain  and  the  system  mentioned  in  Ref.  2  used 


HO 


cathode  followers  as  buffer  amplifiers.   A  way  to  avoid  the 
use  of  such  amplifiers  was  studied  in  this  work. 

1.  Factor  of  Isolation 

A  filter  of  multiple  cells  can  be  constructed  in 
such  a  way  that  the  time  constant  is  maintained  the  same 
through  the  cells  and  the  resistance  is  multiplied  by  a 
factor  of  k  from  cell  to  cell  as  in  Figure  12.   It  can  be 
shown  that  this  is  equivalent  to  a  certain  amount  of  isola- 
tion between  cells  which  depends  on  the  value  of  the  men- 
tioned factor: 

Rp   k  R, 

«=■  =  -5-^  =  k  =  Factor  of  Isolation       (49) 

Rl    Rl 

The  maximum  value  of  k  that  can  be  used  depends  upon 
the  necessary  number  of  cells,  and  the  values  of  the 
capacitors  and  resistors  available. 

2.  Computer  Analysis 

In  the  Appendix  a  program  is  given  which  allows  the 
calculation  of  the  Laguerre  functions  for  any  values  of  a, 
n  and  t.   Figures  13  through  23  show  the  Laguerre  functions 
for  a  =  10   and  for  orders  0  up  to  10. 

The  CORNAP  circuit  analysis  program  was  used  to 
obtain  the  impulse  responses  of  circuits  of  several  number 
of  cells  and  different  values  of  the  factor  of  isolation. 

Impulse  responses  of  circuits  up  to  seven  cells  (n=6) 
were  obtained  because  the  results  had  considerable  errors 
for  circuits  with  a  large  number  of  cells. 
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The  results  showed  that  circuits  with  a  factor  of 
Isolation  of  10  have  exactly  the  same  impulse  response  that 
the  ideal  ones.   The  lowest  factor  of  isolation  that  gave 
reasonable  approximation  to  the  ideal  impulse  responses  was 
found  to  be  5.  Hence,  a  factor  of  isolation  of  10  can  be 
used  whenever  it  is  possible  and  the  use  of  values  less  than 
5  should  be  avoided. 

3.   Laboratory  Results 

Figures  24  and  25  show  pictures  of  impulse  responses 
of  two  kinds  of  circuits:   first  for  a  =  10^  and  a  factor  of 
isolation  of  10  for  orders  0  up  to  3>  and  second  for  the 
same  a  and  k=5  for  orders  0  up  to  5. 

It  can  be  seen  that  the  practical  results  obtained 
agree  with  the  theoretical  considerations  and  computer 
predictions. 

It  was  not  possible  to  build  circuits  with  a  larger 
number  of  cells  because  of  the  characteristics  of  the  avail- 
able pulse  generator,  oscilloscope  and  components. 
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n=3 


Impulse  Responses 
of  Circuits  with 
k=10 


n=2 


FIGURE  24 


n=l 


n=0 
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n=0 


n=l 


n=2 


Impulse  Responses 
of  Circuits  with 
k=5 


n=5 


FIGURE  25         n=4 


n=3 
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VI.   CONCLUSION 

The  realizability  of  the  method  of  Laguerre  type  filters 
or  without  a  "pure"  time  delay  to  obtain  correlation  func- 
tions was  shown  together  with  a  set  of  very  simple  filters 
to  achieve  them  for  low  frequency  signals. 

The  following  things  require  further  investigation: 

Multipliers  and  adder  that  could  be  used  with  the  filters 
presented  in  this  thesis  to  build  a  correlator  as  the  one 
shown  in  Figure  5. 

Factibility  of  building,  through  the  techniques  of  Large 
Scale  Integration,  Laguerre  type  filters  for  several  values 
of  the  base  time  constant  alfa  and  with  a  behavior  close  to 
the  ideal  type  of  filter,  which  means  that  the  circuits  have 
to  use  buffer  amplifiers. 
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APPENDIX 

PROGRAM  TO  COMPUTE  THE  VALUES 
OP  LAGUERRE  FUNCTIONS 


C 

C  MAIN  PROGRAM 

C 

C      THIS  PROGRAM  CALCULATES  THE  LAGUERRE  FUNCTIONS  FOR 

C      ORDERS  0  UP  TO  N  AS  A  FUNCTION  CF  THE  ARGUMENT  X. 

C 

C      THE  PRODUCT  ALFA-^TIME  IS  TAKEN  AS  THE  ARGUMENT. 

C 

C  EASIC    TIME    CONSTANT:KALFA 

C 

C      NUMBER  OF  DISCRETE  VALUES  OF  THE  ARGUMENT :M 

C 

C      INCREMENT  OF  THE  ARGUMENT : DELTA 

C 

DIMENSION  Y(350) 

READ(5tlOO)N 

READ<5, 100JM 

READ(5,101 JDELTA 
C 
C  STARTING    VALUE    OF    ARGUMENT 


REAO(5,101)X 

READC5,  102 JKALFA 

L=N+1 

DO    3C    I=1,M 

CALL    LAP(Y,X,N) 

DAMP=KALFAv£XPl-X) 

DO  2  0  J=1,L 

Y( J) =Y ( J  J  DAMP 

WRITE(6,200i Y( J) 

20  CONTINUE 

X=X+CELTA 

3C  CONTINUE 

100  F0RMATU2) 

101  FORMAT(F5.3) 

102  FORMAT (17) 
200  F0RMATCF9.2) 

STOP 

END 
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c 

C  SUBROUTINE 


C 


SUBROUTINE  LAP 


C 

C         PURPOSE 

C  COMPUTE  THE  VALUES  OF  THE  LAGUERRE  POLYNOMIALS 

C  L(N,X)  FOR  ARGUMENT  VALUE  X  AND  ORCERS  0  UP  TO 

C  N 

C 

C         USAGE 

C  CALL  LAP(Y,X,N) 

C 

C         DESCRIPTION  OF  PARAMETERS 

C  Y    -RESULT  VECTOR  OF  DIMENSION  N+l  CCNTAINING 

C  THE  VALUES  OF  LAGUERRE  POLYNOMIALS  OF 

C  ORDER  0  UP  TO  N  FOR  GIVEN  ARGUMENT  X 

C  VALUES  ARE  ORDERED  FRCM  LOW  TC  HIGH  ORDER 

C  X    -ARGUMENT  CF  LAGUERRE  POLYNOMIAL 

C  N    -ORDER  OF  LAGUERRE  POLYNOMIAL 

C 

c 

SUBROUTINE  LAP(Y,X,N) 
C 

DIMENSICN  Yd) 
C 
C         TEST  OF  ORDER 

Y(l)=l. 

IF(N)1, 1,2 

1  RETURN 
C 

2  Y(2)=l.-X 
IF(N-l)  1,1,3 

C 

C         INITIALIZATION 

3  T=l.+X 
C 

DO  4  1=2, N 

4  Y(I  +  1)  =  Y(  I)-Y(I-1)+Y( I)-(T^Y( I)-Y(  I-11J /FLOAT! I  3 
RETURN 

END 
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